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Abstract 

A pair of coupled erbium doped fiber ring lasers is used to explore the dynamics of coupled 
spatiotemporal systems. The lasers are mutually coupled with a coupling delay less than the 
cavity round-trip time. We study synchronization between the two lasers in the experiment and 
in a delay differential equation model of the system. Because the lasers are internally perturbed 
by spontaneous emission, we include a noise source in the model to obtain stochastic realizations 
of the deterministic equations. Both amplitude synchronization and phase synchronization are 
considered. We use the Hilbert transform to define the phase variable and compute phase synchro- 
nization. We find that synchronization increases with coupling strength in the experiment and the 
model. When the time series from two lasers are time-shifted in either direction by the delay time, 
approximately equal synchronization is frequently observed, so that a clear leader and follower 
cannot be identified. We define an algorithm to determine which laser leads the other when the 
synchronization is sufficiently different with one direction of time shift, and statistics of switches in 
leader and follower are studied. The frequency of switching between leader and follower increases 
with coupling strength, as might be expected since the lasers mutually influence each other more 
effectively with stronger coupling. 
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The main goal of this paper is to explore the synchronization of mutually 
delay-coupled spatiotemporal systems and develop techniques to identify the 
evolving phase relationships between them. We illustrate our ideas in the specific 
case of mutually coupled fiber ring lasers, which exhibit a remarkable ability to 
adapt their dynamics so that much of the time there is no clear leader or follower, 
even though signals propagate from one system to the other with an appreciable 
time delay. The lasers are described by a system of stochastic delay differential 
equations, so as to include the effect of spontaneous emission in the erbium doped 
fiber amplifiers that constitute the active media. We use statistical measures to 
study phase and amplitude synchrony in model and experimental delay-coupled 
systems from their time series. 

One fascinating area of dynamics in nature is the study of how systems respond to each 
other due to their interactions. When systems behave dynamically, interactions between 
them may cause the systems to operate in a similar, or coherent, manner. If the measured 
signals are similar enough over time, the correlated motion is termed synchronized. There 
now exist several excellent reviews on synchronization dynamics in the literature, such as 
Refs. , and they cover a wide range of applications from many fields of science. 

In general, synchronization between interacting systems may be quantified by examining 
and comparing the output time series from each dynamical system. Although the dynamics 
may be considered for the general coupling between N systems, in this paper we restrict the 
number of dynamical systems to iV = 2. Then we may consider only two coupling schemes; 
either mutual (bidirectional), or unidirectional. 

Suppose Xj : R — > R n , i = 1, 2 denote the vector output time series measured from each 
dynamical system. Then several types of synchronization may be classified depending on the 
type of coherence measure. The systems are in complete synchronization if Xi(i) = X 2 (t). 
Complete synchronization occurs in coupled phase oscillators as well as in coupled chaotic 
oscillators pi, In this case, amplitudes and phases are identical. 

If the amplitudes are uncorrelated but the phases are locked, or entrained, between the 
two signals, then the systems are said to be in phase synchrony [7J. One other type of syn- 
chronization deals solely with the unidirectional coupling between two oscillators of drive 
and response type, and is termed generalized synchronization 0. In generalized synchro- 
nization, there exists a functional relationship between the drive and response, where there 
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might exist a function F such that Xi(t) = F(X2{t)). In a more general setting, this may 
also be thought of as a generalized entrainment in dynamics, whereby one system is en- 
trained functionally to another. Many examples of entrained systems occur in singularly 
perturbed problems, and specifically in systems with highly different time scales jjj. 

If there is a parameter mismatch or noise in the dynamical systems, complete synchroniza- 
tion may not be possible, and other measures of synchronization are needed. One possible 
example of phase synchronization occurs when the amplitudes are correlated but locked in 
phase at a value other than zero. Such a system exhibits lag synchronization when for some 
r > we have Xi(i) = X 2 (t — t). That is, the two outputs of the coupled dynamical systems 
appear shifted in time with respect to one another. For mutually coupled chaotic systems 
such as Rossler attractors that are mismatched in frequency, lag synchronization is one of 
the routes to complete synchrony as coupling is increased (l0|. It should be noted that lag 
synchronization may occur without the presence of delay in the coupling terms. 

On the other hand, if delay is introduced into the coupling terms to model finite time 
signal propagation, then synchronized behavior may still occur. When there is a clear time 
lag between the delay-coupled dynamics, the systems are said to exhibit achronal synchro- 
nization. Achronal synchronization exhibits a clear leading time series which is followed by 



a lagging time series. Heil et al. 11] showed the existence of achronal synchronization in 
a delay-coupled semiconductor laser experiment, as well as in a single mode model of the 
delay-coupled lasers in which stochastic effects modeling spontaneous emission are included. 
The time series shows a clear leader with delay equal to the coupling delay time. Other 
groups have also considered leader-follower synchronization in single-mode semiconductor 
equations jl^ . 

One of the non-intuitive facets of interacting systems which synchronize with delay is 
that of anticipation in systems with short time delay. First observed in unidirectionally 
coupled systems, in contrast to lag synchronization, anticipatory synchronization occurs 
when a response in a system's state is replicated not simultaneously but anticipated by the 
response system jlH Q, Q|. An example of anticipation in synchronization is found in 
coupled semiconductor lasers [16]. Here, the author followed Ref. 14] in the design of a 



unidirectional coupling arrangement, in which two single mode semiconductor lasers with 
delayed optical feedback and delayed injection coupling were modeled. Cross-correlation 
statistics between the two intensities showed clear maxima at delay times consisting of 
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the difference between the feedback and the coupling delay. Anticipatory responses in the 
presence of stochastic drives, equally applied to transmitter and receiver, have been observed 



in models of excitable media as well 



12|. 



Given that both lag and anticipatory dynamics may be observed in delay-coupled systems 
which are deterministic, it is natural to ask whether the systems may exhibit coexisting 
features. In single mode mutually coupled lasers, this is indeed the case In fact, in 

the absence of noise (the authors include a noise term in their model, but turn it off for 
the simulations), switching between leading and following state is observed. In a similar 
model with spontaneous emission included as a noise source, theory and experiment have 
exhibited achronal synchronization with switching between leader and follower. The 
authors conjecture that for semiconductor lasers, changes in leader-follower roles may occur 
during sharp dropouts of the laser intensity. 

In the above discussion, it is clear that the role of signal propagation time is important 
in the dynamics of leader and follower in coupled systems. When noise is sufficiently large, 
it is not easy to detect changes in leader and follower statistically. This is espec ially the 
case when the systems are stochastic and close to a synchronized state. In Ref. the 
authors compute statistics on the phase differences and show that they might correspond 
to a random walk. In this paper, we develop methods to extract statistical behavior of the 
switches in leader and follower in mutually coupled fiber lasers. 

The chaotic dynamics of fiber ring lasers have been studied in the past. An experiment 
on the coupling between polarization modes was set up and modeled using delay differen- 
tial equations i n | 2C |. Other experiments on synchronization with fiber lasers have been 
reported in j^li I22I 23L and noise- induced generalized synchronization in fiber ring lasers 
has been reported in [24j. Modeling the ring laser yields a system of equations which con- 
sists of coupled difference and differential delay equations. To obtain better agreement with 
experiment, it was found that inclusion of spontaneous emission effects was necessary in the 
modeling, which resulted in a stochastic difference-differential system of equations [2j| , and 
it is this approach we follow here. 

For mutual coupling of two ring lasers with delay, we analyze both experimental results 
and a delay differential equation model of the system. Synchronization with delay occurs 
and increases with coupling strength. Generally, approximately equal synchronization is 
observed between the two lasers when time-shifted in either direction by the delay time, so 



5 



that a clear leader and follower cannot be identified. We define an algorithm to determine 
which laser leads the other when the synchronization is sufficiently different, and statistics 
of switches in leader and follower are studied. 

We introduce the experimental setup and corresponding model in Sec-HHiTl The dynamics 
of the system are outlined in Sec. 11111 Synchronization of the coupled lasers is discussed in 
Sec. |IV] and switching of leader and follower in Sec. |V| 

I. EXPERIMENTAL METHODS 

The experimental setup consists of two erbium doped fiber ring lasers (EDFRLs) coupled 
with two passive coupling lines. Each EDFRL has approximately f 7 m of erbium-doped 
fiber, the active medium, and approximately 29 m of passive single (transverse) mode fiber, 
making the total length of each cavity approximately 46 m. The lengths of the cavities are 
matched within 1 cm of each other. Though the doping density and the lengths of the active 
media are the same in both lasers, defects and imperfections in the fiber make the lasers 
nonidentical. The erbium ions in the active medium are pumped with identical 980 nm 
semiconductor lasers at a pump power of 120 mW. The lasing threshold for both lasers 
is approximately 20 mW. This results in approximately 1 mW of light circulating within 
each ring. An optical isolator is inserted within each ring cavity to ensure unidirectional 
propagation within the rings. 

Each laser contains 4 fiber-optic evanescent field couplers placed at similar locations in 
both rings. These consist of two 70/30 couplers which input and output light between the 
lasers, one 90/10 coupler for monitoring, and a 95/5 coupler as an extra port for applications 
that are beyond the scope of this paper. The locations of the couplers are shown in Figure 
^ The ports of the couplers not in use are angle cleaved to ensure that there are no back 
reflections and were monitored to ensure that light was propagating in the correct direction 
within the cavities. 

The lasers are connected via two injection lines, which consist of passive single mode 
optical fiber, one splitter, and a variable attenuator. In this configuration, we have the ability 
to monitor and control the injection amplitude between the lasers, through the splitter and 
variable attenuator respectively, and to observe it on an oscilloscope. The coupling strength 
is defined to be the ratio of the power of light in the injection line to the power of light 
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FIG. 1: Experimental setup for the coupled fiber laser system. 

in the source ring. (Note that both lasers have been adjusted to have the same power.) 
The lowest coupling strength that can be resolved in our system is 0.0001. The injection 
lines are approximately 9 m long, corresponding to a travel time between the two lasers of 
approximately 45 ns, and again they are matched within 1 cm. In the experiments presented 
here, the two coupling strengths are always the same (symmetric coupling), though the 
electric field from each laser undergoes different phase and polarization changes due to fiber 
imperfections along their separate paths. 

The electric field intensity of each laser is monitored using a 125 MHz bandwidth pho- 
todetector and a 1 GSample/s digital oscilloscope. The optical spectrum of each laser is also 
monitored, and the spectra of the uncoupled lasers are matched to within 1 nm. 



II. DESCRIPTION OF MODEL 



We use the model of |25j] with the modifications introduced in [2y, |22[ for coupled fiber 
ring lasers. Only a single polarization mode in each laser is considered here. The model is 
characterized by the total population inversion W(t) (averaged over the length of the fiber 
amplifier) and the electric field E(t) in each laser. Two delay times occur in the model: tr, 
the cavity round-trip time, and tj, the delay in the coupling between lasers. The equations 



for the model dynamics are as follows: 




(1) 



(2) 



where 



E[^(t) = E 1>2 {t - t r ) + nE 2A (t - r d ). 



(3) 



Ej(t) is the complex envelope of the electric field in laser j, measured at a given reference 



laser j and optical coupling with the other laser. Time is dimensionless, measured in units 
of the decay time of the atomic transition, 711. The active medium is characterized by 
the dimensionless detuning otj between the transition and lasing frequencies and by the 



and Nq the population inversion at transparency The ring cavity is characterized by its 
return coefficient R, which represents the fraction of light remaining in the cavity after one 
round-trip, and the average phase change A<ft = 27rnL p /\ due to propagation of light with 
wavelength A along the passive fiber of length L p and index of refraction n. Energy input is 
given by the pump parameter q. The electric field is perturbed by a complex Gaussian noise 
source £j with standard deviation D. Coupling between the lasers is characterized by the 
coupling strength k, which was varied from to 0.014. The same k is used in both lasers 
(symmetric mutual coupling). Values of the parameters in the model, similar to those in the 
experimental system, are given in Table |1J 

Eqns. mm consist of a delay differential equation for W(t) coupled to a map for E(t). 
We integrated Eqn. El numerically using Heun's method. The time step for integration was 
tr/N, where N = 600. This step size corresponds to dividing the ring cavity into N spatial 
elements. 

Because of the feedback term Ej dh (t) in Eqn. ^ one can think of Eqn. ^as mapping the 
electric field on the time interval [t — Tr, t] to the time interval [t, t + tr] in the absence of 
coupling (k = 0). Equivalently, because the light is traveling around the cavity, Eqn.m ma ps 
the electric field at all points in the ring at time t to the electric field at all points in the 



point inside the cavity. E® 



(t) is a feedback term that includes optical feedback within 



dimensionless gain T = ^aL a N 0l where a is the material gain, L, 



a 



the active fiber length, 



S 



TABLE I: Parameters used in the coupled fiber laser model. 



Parameter 


Value 


Units 


Description 


R 


0.4 




output coupler return coefficient 


a 


2.03 x 10~ 23 


m 2 


material gain coefficient 


L a 


15 


m 


length of active fiber 


L p 


27 


m 


length of passive fiber 


No 


1Q 20 


m -3 


transparency inversion 


r 


0.0152 




dimensionless gain 


a x 


0.0352 




detuning factor, laser 1 


a 2 


0.0202 




detuning factor, laser 2 


n 


1.44 




index of refraction 


A 


1.55 x 10~ 6 


m 


wavelength 


A0 


1.58 x 10 8 




average phase change 


D 


0.02 




standard deviation of noise 


q 


100 




pump parameter 


7|| 


100 


s- 1 


population decay rate 


tr 


201.6 x 10~ 9 


s 


cavity round-trip time 


Td 


45 x 10~ 9 


s 


delay time between lasers 


K 


0-0.014 




coupling strength 



ring at time t + tr. We can thus construct spatiotemporal plots for E(t) or the intensity 
I(t) = \E(t)\ 2 by unwrapping E(t) into segments of length tr. 

The experimental measurements were of the intensity of light from each laser after passing 
through a 125 MHz bandwidth photodetector. To correspond with the experiment, we com- 
puted intensities from model and applied a low pass filter with f = 125 MHz, multiplying 
the Fourier transform by the transfer function 

2 



G 




7 

Jo, 



f 

Jo 



(4) 



It should be noted that the time step for integration is an order of magnitude smaller than 
the smallest timescale allowed by the filtering function. We expect that the results are 
insensitive to the integration time step as long as it is sufficiently small. Test runs with a 
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finer mesh, iV = 1200 time steps per round-trip, yielded similar behavior to that described 
here. 

All subsequent analysis of the model properties is based on the filtered intensity. 

III. EFFECT OF COUPLING ON SYSTEM DYNAMICS 

In this section we discuss the types of dynamics observed in simulations and the exper- 
imental system. After transients have died out, each laser settles into a pattern that is 
approximately repeated every round-trip. The pattern may shift or change over time inter- 
vals of tens or hundreds of round-trips. Due to the presence of noise in the model, even the 
uncoupled lasers (re = 0) are not fully periodic for these parameter values. 

Typical spatiotemporal plots and time traces for the model are shown in Figure 121 Spa- 
tiotemporal plots are constructed by displaying each round-trip as a row in the diagram, 
colored according to the laser intensity, with subsequent round-trips forming subsequent 
rows. 

The system displays approximately steady behavior for small coupling (k < 0.003) 
(Fig. |2K)- As the coupling increases (k > 0.003), it then enters a regime where pulsing 
occurs, in which the intensity alternates between pulsing and dropping to a low value while 
the inversion builds up (Fig. Eb)- Although the spatiotemporal plot may appear more peri- 
odic during the pulses, the intensity is actually highly irregular. When the coupling increases 
further (k > 0.007), the system leaves the pulsing regime and displays complicated behavior. 
Traveling wave solutions are commonly observed (e.g., Fig. |2b). As the coupling increases, 
the traveling waves become less prominent (e.g., Fig. EJi). Time traces of the intensity 
are also shown (Figure Et-h) for two adjacent round-trips. Approximate repetition of the 
intensity pattern from one round-trip to the next can be seen. 

Spatiotemporal plots for the experiment are shown in Figure 01 Because the round-trip 
time is not an integer multiple of the 1 ns sampling time, spacetime plots were constructed 
by the following procedure. The experimental intensity time trace was expanded to 10 times 
as many time points by linear interpolation. Subsequent round-trips were aligned by shifting 
the relative position of the rows to maximize correlation between the rows. Spatiotemporal 
plots for the experiment are normalized so that the intensities range from to 1. 

The experimental system displays similar behavior to the model. When k = 0, approx- 
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FIG. 2: (a)-(d): Spatiotemporal plots of Laser 1 intensity for the model, [(a) n = 0, (b) k = 0.005, 
(c) k = 0.009, (d) k = 0.011.] (e)-(h): Model-'- intensity vs. time for two round-trips. [Same k 
values as in (a)-(d). Note that (f) shows typical behavior in the low intensity regions of the pulsing 



imately steady behavior is observed (Figure Ek)- When very small coupling is turned on, 
higher frequency structures emerge in the spatiotemporal pattern (Figure Eb)- Complicated 
spatiotemporal patterns are seen at stronger coupling (Figure Et)- Pulsing was not observed 
in this experimental study. However, pulsing in other fiber ring laser systems has been com- 
monly reported in the literature (e.g., [28]). As in the model, the experimental time traces 
show the approximate repetition of the intensity pattern in subsequent round-trips (Figure 

BH-f). 



IV. SYNCHRONIZATION 



We next consider synchronization between the two lasers. Because the lasers are mutually 
delay-coupled, with a delay time r d , we expect the coupling to cause lag synchronization in 
the lasers so that Ii(t) = l2(t — r d ) or I 2 (t) = Ii{t — T d ). In FigurelH we display time traces of 
the two lasers with an offset of r d . For both the experiment and the model, synchronization 
can be seen with a time shift in either direction. 

We studied the synchronization quantitatively using two metrics, the cross- correlation C 
for amplitude synchronization and the mean phase coherence R for phase synchronization. 
Let Cyi denote the correlation computed for Laser 1 leading Laser 2 and C21 the correlation 
for Laser 2 leading Laser 1. Thus for time series of length K time points, 



a 



ab 



1 1 

a a a b K 



K 



^[Uh ~ T d ) - {I^Wj) - {I b % 



(5) 



0=1 



where (•) denotes averaging and a a , o~b are standard deviations in intensity for Lasers a and 
b. 

We define Ru and R21 similarly. The phase of a time series is defined as the Hilbert phase 



<p(t), computed from the Hilbert transform of the intensitv 



two time series of length K is then defined as follows 

K 



29 



30 



he mean phase coherence for 



R, 



ab 



K 



A<t>, 



ab,j 



(6) 



where 



A0a6,j = <fia(tj - T d ) - <f) h (tj) (7) 

is the phase difference between the jth points of the time series when Laser a leads Laser b. 
R ranges between (no synchronization) and 1 (complete phase synchronization). 
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FIG. 3: (a)-(c): Spatiotemporal plots of Laser 1 intensity for experiment (intensities rescaled to 
have equal ranges), [(a) n = 0, (b) k = 0.000285, (c) k = 0.0228.] (d)-(f): Experimental intensity 
vs. time for two round-trips. [Same n values as in (a)-(c).] 
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FIG. 4: Time traces of intensity (arbitrary units). In each plot, the center bold curve is intensity 
vs. time for a single round-trip in Laser 1. The lighter curves above and below it are intensity 
vs. time for the same time interval in Laser 2, but plotted with the curve shifted forwards or 
backwards in time by the delay time r^. The correspondence between the bold and light curves 
indicates synchrony between the two lasers with a shift of r^. (a) Experimental results for n = 
0.0114, (b) model results for « = 0.011. 

The synchronization metrics were calculated for each round-trip, with each laser leading. 
We divided the time series into round-trips and shifted each short time series of length tr by 
the delay in either direction. (Hanging ends were omitted, so that only tr — t& was used 
in the calculation.) C and R were then computed. When analyzing the experimental data, 
we compensated for the discrete sampling by using linear interpolation to expand the data 
to 10 times as many time points and then calculating C and R for a 2.5 ns range of coupling 
delays (centered around the estimated r d ). The maximal Ci 2 , C 2 i, -R12, and _R 2 i values for 
each round-trip (maximized over possible coupling delays) were recorded. Sample results 
for phase synchrony are shown in Figure El The synchronization fluctuates over time, and 
there are changes in which laser leading gives the maximum synchronization. The amplitude 
synchronization behaves similarly. 

To determine the mean amplitude synchronization (C) for a given coupling value, we 
take the maximum of C12 and C21 for each round-trip and then average over all round- 
trips. We similarly define the mean phase synchronization (R) over all round-trips. Mean 
synchronizations were computed for each value of the coupling constant k. Experimental 
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FIG. 5: Phase synchrony vs. round-trip number for experiment with k = 0.0114 coupling. Blue 
curve: synchrony computed with Laser 2 leading (i?2i); red curve: synchrony computed with Laser 
1 leading (R^)- Locations of switches identified by our algorithm are indicated by arrows. 

data were averaged over 5-10 separate runs of 4500 round-trips each. Simulations were 
averaged over 6 runs of 4 x 10 4 round-trips each (with runs separated by at least 0.1 s of 
simulated time to insure independence). Results are shown in Figure E3 Error bars are the 
standard deviations in C and R over all round-trips, so they are a measure of fluctuations 
in synchrony. 

In both the model and experiments, amplitude and phase synchronization increase as 
the coupling increases. However, the shapes of the synchronization vs. coupling curves 
are different. The experimental system is fairly well synchronized even at the smallest 
nonzero coupling, while the model synchronizes more gradually. The model also exhibits 
more fluctuations in synchrony, as exemplified by the larger standard deviations. For both 
model and experiment, the amplitude synchronization (cross-correlation) is typically larger 
than the phase synchronization. 

V. SWITCHING 

We next consider which laser leads the other during synchronization, and switches in 
the leader /follower. In Figure the two curves are the synchrony computed with Laser 1 
leading (Ru) and with Laser 2 leading (i?2i)- Examining the figure qualitatively, we see 
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FIG. 6: Average correlation between lasers vs. coupling k for experiment (a) and model (b) and 
average phase synchrony between lasers vs. k for experiment (c) and model (d). Error bars indicate 
standard deviation over round-trips. 

that at certain times the synchrony is clearly greater with one laser leading than with the 
other (e.g., round-trips 1300-1900), and at other times it is more difficult to determine which 
laser is leading (e.g., round-trips 2500-2900) because the difference in synchrony with time 
shifts in either direction is small compared to the fluctuations in the synchrony. The latter 
portions of Figure El are enlarged in Figure We develop an algorithm to compute the 
leader and follower. 

For each round-trip, we compute the synchronization with either laser leading. Let ACj = 
Ci2{i) — C2i{i) be the difference between the correlations for each round-trip % within a run, 
and let gq be the standard deviation of the ACj. When |ACj| > hoc for some cutoff factor 
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FIG. 7: Phase synchrony vs. round-trip number for experiment with 0.0114 coupling. Blue curve: 
synchrony computed with Laser 2 leading (R21); re d curve: synchrony computed with Laser 1 
leading (R\2)- Locations of switches identified by our algorithm are indicated by arrows. 

h, we say that the synchronization is substantially greater with one laser leading than with 
the other, and a clear leader and follower can be identified. Leader and follower for phase 
synchronization are defined similarly. Figure |H1 shows the fraction of round-trips for which 
a leader and follower can be identified for a range of cutoff factors h. Similar curves are 
obtained for the model and experiment for most values of k. We select h = 1 to use as the 
cutoff factor in this study. For most values of k, for both the model and experimental data, 
the choice h = 1 leads to approximately 30% of round-trips having an identified leader and 
follower for amplitude and phase synchronization. 

After leader and follower have been computed, we next search for switches. The following 
algorithm is used to identify changes in leader and follower. Let {n^} be the set of round- 
trip numbers for which |AC n J > hac- Whenever the sign of AC ni is opposite the sign of 
AC ni _j, a switch is identified at round-trip rii. 

To obtain statistics over several runs, let be the total number of round-trips and N sw ^c 
total number of switches in amplitude synchronization observed. Then the average number 
of round-trips between switches is T c = N t /N sw ^ C - We assume square root error in the 
number of observed switches, leading to an error of N t /N^ c in the number of round-trips 
between switches. Switches in phase synchronization are located by the same procedure. An 
example of switching is given in Figure El and Figure [3 locations of switches identified by 
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FIG. 8: Fraction of round-trips for which leader and follower can be identified vs. number of 
standard deviations used in cutoff (h) for model with k = 0.01. 

our algorithm are marked by arrows. Note that switches are defined only when the leader 
and follower change; for example, no switch is identified around round-trip 2000 when R\2 
and i?2i come together with no defined leader and then separate with the same leader and 
follower as before. 

Switches were studied in the model and experiment for a range of coupling strengths, 
using the same data sets as in Section IIVI Results are shown in Figure M We find that 
the time between switches in phase synchronization decreases as the coupling increases for 
both the experiment and model. Results for the amplitude synchronization are less clear. 
For the model, switching in amplitude synchronization follows the same trend as for phase 
synchronization. However, the experimental system exhibits fewer switches in amplitude 
synchronization at larger coupling strengths. 

VI. DISCUSSION AND CONCLUSIONS 

We compared experimental measurements with model predictions for a pair of mutually 
coupled fiber ring lasers as the coupling strength was varied. The coupling has similar effect 
on the dynamics in both model and experiment. Quantitative measures of amplitude and 
phase synchronization were made, the latter by using the Hilbert phase. Approximately 
equal synchronization was found for forward and reverse time shifts equal to the coupling 
delay. We notice that the lasers spend considerable periods of time in synchronized states 
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FIG. 9: Average number of round-trips between switches in leader/follower as determined by 
cross-correlation [(a) experiment, (b) model] and by phase synchrony [(c) experiment, (d) model]. 

with no discernable leader or follower, since the synchronization values are very similar for 
time shifts in either direction. This is in contrast to the mutually coupled system with 
delay but no optical feedback described in 19], where leader and follower could be identified 
when the cross-correlation was computed for a short time series. We defined an algorithm 
to locate subtle switches of leader /follower between the two lasers. Our statistical measures 
allow a leader and follower to be identified when the synchronization is sufficiently larger in 
one direction of the time shift. 

Although the trend of increasing synchronization with increasing coupling was consistent 
between model and experiment, the shape of the synchronization vs. coupling curves differed, 
with the experimental system synchronizing more quickly. In addition, the fluctuations in 
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synchronization were generally larger in the model. The model we have used includes only 
one polarization mode in each laser. Future extensions to the model should include both po- 
larization modes. This may explain the discrepancy between the experimental observations 
and computational prediction in the shapes of Figure El 

In |l9[, the authors note that for long time series, forward and reverse time shifts 
yielded similar cross-correlations, although leader and follower could be identified from cross- 
correlations of short time series. Our cross-correlations were computed over the natural time 
scale of a round-trip, and similar results are obtained for shorter time series with length equal 
to the delay time (data not shown). Our observation that the system often does not have a 
clear leader and follower seems robust for reasonable choices of window size for averaging, 
although further study of appropriate window size may be conducted. 

The method for choosing the cutoff to identify clear regions of leader and follower is 
based on an internal definition for a given time series. The choice of the cutoff factor h to 
be one standard deviation is somewhat arbitrary, but preliminary results indicate that the 
trends of the switching frequency vs. coupling are preserved for different cutoff factors. The 
algorithm for identifying switches may in the future be applied to other systems in which 
the difference in synchronization for forward and reverse time shifts is more substantial. 

We typically observe a higher frequency of switching at higher coupling values. There 
are very similar trends in the frequency of switches of leader /follower for the phase synchro- 
nization in experiment and theory. This effect might be expected since the lasers mutually 
influence each other more effectively with stronger coupling and thus can respond to each 
other more quickly. However, we cannot explain the differences in corresponding plots for 
the amplitude synchronization. Further study of this phenomenon is needed. 

In this study, we considered variations in the coupling strengths only. The noise strength 
for the numerical computations was chosen somewhat arbitrarily. The main guiding prin- 
ciple was to have a reasonable correspondence between the uncoupled lasers in theory and 
experiment. Future studies will consider changes in the noise strength. Another parameter 
to be varied is the detuning between lasers, because its value is not known exactly for the 
experimental system. The inclusion of a second polarization mode will introduce another 
set of parameters whose variation may also be studied. 

We are indebted to Jordi Garcia Ojalvo for developing the original version of the model 
and the simulation code, and we acknowledge David DeShazer's assistance in generating 



20 
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Naval Research. LBS is currently a National Research Council post doctoral fellow. EAR 
was supported in part by a National Science Foundation fellowship. 
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